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The polarization ol the cosmic microwave background (CMB) is widely recognized as a potential 
source of information about primordial gravitational waves. The gravitational wave contribution 
can be separated from the dominant CMB polarization created by density perturbations at the 
times of recombination and reionization because it generates both E and B polarization modes, 
whereas the density perturbations create only E polarization. The limits of our ability to measure 
gravitational waves are thus determined by statistical and systematic errors from CMB experiments, 
foregrounds, and nonlinear evolution effects such as gravitational lensing of the CMB. Usually it 
is assumed that most foregrounds can be removed because of their frequency dependence, however 
Thomson scattering of the CMB quadrupole by electrons in the Galaxy or nearby structures shares 
the blackbody frequency dependence of the CMB. If the optical depth from these nearby electrons 
is anisotropic, the polarization generated can include B modes even if no tensor perturbations are 
present. We estimate this effect for the Galactic disk and nearby extragalactic structures, and 
find that it contributes to the B polarization at the level of ~ (1-2) x 10 _4 ^K per logarithmic 
interval in multipole t for £ < 30. This is well below the detectability level even for a future CMB 
polarization satellite and hence is negligible. Depending on its structure and extent, the Galactic 
corona may be a source of B-modes comparable to the residual large-scale lensing B-mode after the 
latter has been cleaned using lensing reconstruction techniques. For an extremely ambitious post- 
Planck CMB experiment, Thomson scattering in the Galactic corona is thus a potential contaminant 
of the gravitational wave signal; conversely, if the other foregrounds can be cleaned out, such an 
experiment might be able to constrain models of the corona. 

PACS numbers: 98.80.Es, 95.30.Gv, 98.35.Gi 



I. INTRODUCTION 

Recent observations of the cosmic microwave back- 
ground (CMB) anisotropies have confirmed several of the 
predictions of the simplest inflationary models Q, 0, • 
They have (in combination with other datasets) shown 
that the universe is close to spatially flat, with Sl tot = 
1.02 ± 0.02 0; that the perturbations in the CMB 
are close to Gaussian 0; that the temperature and E- 
polarization power spectra match those expected from 
adiabatic initial conditions 0; and that the perturba- 
tions are close to scale-invariant with n s « 1 and a s w 
. Thus there is now interest in testing the last of the 
major predictions of inflation: the existence of a roughly 
scale-invariant background of tensor (gravitational wave) 
perturbations. The amplitude of the tensor perturba- 
tions is typically described with the ratio of the temper- 
ature quadrupoles in the CMB T/S = Cf nsor / 'C^ aI '. 
The value of T/S depends on the specific inflationary 
model through the relation T/S = 14/(3.7 x 10 16 GeV) 4 , 
where V* is the energy density at the time during infla- 
tion when observable scales in the CMB (10 -4 < k < 
10~ 2 h Mpc^ 1 ) exited the horizon. A detection of the 
tensor perturbations would bolster confidence in infla- 
tion, as well as providing a measurement of V* . 

If T/S is large, then the tensor perturbations add to 



the I < 100 temperature power spectrum of the CMB, 
and the nondetection of this feature in the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) data sets an upper 
limit of T/S < 0.7 0- However the use of the temper- 
ature power spectrum to constrain T/S suffers from the 
large cosmic variance error bars on the low CMB mul- 
tipoles, which impose a fundamental limit on measure- 
ments of T/S from Cj T . In addition, these measure- 
ments are model-dependent in the sense that modifica- 
tions to the scalar power spectrum can produce excess 
power at £ < 100 without invoking tensors (for example, 
T/S is partially degenerate with the running of the spec- 
tral index ot s in the current data 0). An alternative is to 
use the CMB polarization, which can be decomposed into 
E and B modes 0, . Since the scalar density fluctua- 
tions of multipole £ have parity (— 1) , whereas the E and 
B polarization modes have parity (— I) 1 and — (— 1)^ re- 
spectively, it follows that the scalars can produce only E 
polarization in linear perturbation theory. However ten- 
sor perturbations consist of both positive and negative 
parity modes for every value of £ > 2, so the tensors can 
contribute to both E and B. The measurement of T/S 
through £?-mode polarization thus is not limited by the 
cosmic variance associated with the scalars, and is not de- 
generate with any features in the scalar power spectrum; 
instead it is limited by the observations, foregrounds, and 
nonlinear processes that can create i?-mode polarization 
from scalar initial conditions. 
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The dominant nonlinear process that generates B- 
mode polarization in the CMB is gravitational lensing 
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of the E polarization generated during recombination. 
On the large angular scales relevant for tensor B-mode 
searches, the lensing B-mode appears as white noise with 
Cf B w 2.4 x lCrVK 2 ; this would set a detectability 
limit of (T/5') m i n = 6 x 10~ 5 for a full-sky experiment 
if t = 0.17 and (T/S) mi „ = 2x 10~ 5 if r = 0.07. (The 
optical depth to reionization r is relevant since this de- 
termines the amplitude of large-scale polarization.) How- 
ever, lensing results in higher-order correlations between 
the E and B polarization modes that allow one to "clean" 
out the lensing 5-mode [H El [ljj. The amount by 
which one can clean the B-mode using the iterative tech- 
nique of Ref. Q was computed by Ref. as a function 
of the detector noise and beam size; for the most sen- 
sitive experiment considered (0.25 /xK arcmin noise and 
2 arcmin beam) , the B-mode polarization can be cleaned 
down to Cf B w 5.8 x 10~VK 2 , which reduces (T/S) min 
by a factor of 40. 

The polarized foregrounds contain both Galactic and 
extragalactic components. The extragalactic foregrounds 
such as point sources are approximately white noise and 
hence are largest on small angular scales, whereas the 
Galactic synchrotron and dust emission are expected to 
have a significant polarization at low t. These can be at 
least partially cleaned using their frequency dependence, 
which differs from the blackbody signature expected from 
tensors. However it is also possible to generate polariza- 
tion via Thomson scattering of the CMB quadrupole by 
free electrons in the Galaxy or other nearby structures. 
While this polarization signal is extremely small and is 
well below the predicted level of the synchrotron and/or 
dust foregrounds at all frequencies, it cannot be cleaned 
using frequency information since it is a blackbody sig- 
nal. It can contain both E and B polarization if the 
distribution of scattering electrons is anisotropic, hence 
it is a potential foreground for tensor B-mode searches. 
The primary purpose of this paper is to investigate this 
Thomson polarization from the local universe. 



II. GALACTIC CONTRIBUTION TO 
POLARIZATION 

A. Basic model 

The polarization Stokes parameters P = (Q,U) pro- 
duced by Thomson scattering are most easily determined 
in the line-of-sight formalism |l6j | as 



and 



/7* />oo ^ 

P(n)= / g(r,h) £ [T 2m (rn) - V6E 2m (rh)} 



xYf m (n)dr, 



where 



T 2m (rn) = J Y* m (h')T(n')\ rh d 2 h' 



(1) 



(2) 



E 2m (rn) = J Yf„I(n')P(n')|,nrf 2 n' 



(3) 



are the temperature and polarization quadrupole mo- 
ments of the CMB at position rn relative to the observer, 
and 



l(r, n) 



dr(r, n) 
dr 



(4) 



is the Thomson visibility function, which depends on the 
optical depth r(r, n) to distance r; we have trivially gen- 
eralized from Ref. [l|| to allow r to depend on the direc- 
tion n. Here Yi m are the spherical harmonics and Yf m 
are the tensor spherical harmonics. (We normalize these 
to J YfJ[Yf m (n) d 2 n = 1 over the whole sphere, with 

the inner product P|P 2 = QIQ2 + U?U?:, this is con- 
sistent with Ref. and with CMBFast 17] but differs 
by a factor of \[2 from Ref. 10].) If we break the line- 
of-sight integral into a piece contained within our galaxy 
(r < R = 100 kpc) and a piece external to the galaxy 
(r > R), we find to first order in t(R, n): 

P(n) = [1 - t{R, n)]P cosmic (n) + P ga i(n), (5) 

where P cosm i C (n) is the polarization that we compute if 
we neglect the Galactic Thomson scattering and 



P gal (h) = -^r(R 7 h) J2 ^ 

m= — 2 



V6£ 2m )Yf m (n), 



(6) 

(We have assumed that the cosmic quadrupoles 
T 2TO (rn) — y/&E2m{f"n) at position rn inside the Galaxy 
can be replaced by their values at the observer, T 2m — 
V6E2m, which is a good approximation since the cos- 
mic quadrupole is dominated by perturbation modes 
with wavenumber k <C R .) The optical depth is 
related to the electron distribution through the result 
dr I dr = n e aT/(l + z) where n e is the electron density 
and <7t = 6.65 x 10 - 25 cm 2 is the Thomson cross section. 

If t(R, n) were independent of n, i.e. if the Galactic 
optical depth were the same in all directions, then the 
only effect of the Thomson scattering would be to (i) 
suppress the CMB power spectrum because of the factor 
of 1 — t(R, n), and (ii) generate a small amount of £ = 2 
-E-mode through P ga ;- However if t(R, n) is anisotropic 
with fluctuations St/t ~ 1, which is the case in the real 
Galaxy we inhabit, then P ga i becomes a mixture of E 
and B modes that is generically of order r times the cos- 
mic temperature+polarization quadrupole. In the rest 
of this paper, we will neglect -E 2m compared to T 2m in 
Eq. ©, since even for the WMAP r = 0.17 cosmology 
we have C EE = 0.08 ^K 2 , whereas the observed temper- 
ature quadrupole in the map of Tegmark et al. [18j is 
C 2 TT = 200/iK 2 . 

Since this temperature quadrupole has already been 
measured, we can compute the induced polarization if 
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we have a map of the Galactic electron distribution. In 
this paper we compute the polarization P ga i using the 
model of Cordes & Lazio for the Galactic distribution 
of free electrons. This model is based on pulsar disper- 
sion (DM) and scattering (SM) measures as well as SMs 
of other Galactic and extragalactic radio sources, and is 
appropriate for the electron content of the Galactic disk; 
the predicted optical depth in the direction of the Galac- 
tic poles is 4 x 1(T 5 (North) or 5 x 1(T 5 (South). This is 
somewhat larger than the optical depth of ~ 3.4 x 10~ 5 
predicted by the older Taylor & Cordes model po|. 

The input data to the Cordes & Lazio model is sparsely 
sampled, with a typical density of pulsars at high Galac- 
tic latitude (|6| > 20 deg) of n ~ 6sr _1 , and small-scale 
structure in the electron density map is thus expected to 
be lost |2]J. As a crude estimate, we note that the num- 
ber of modes per steradian is ~ ^^ ax /47r; setting this 
equal to n implies that structures smaller than € max ~ 9 
should be lost. (The 14 non-pulsar SMs at high Galactic 
latitude do not significantly alter this conclusion.) The 
effective ^ max at high Galactic latitude may be even less 
than this because some of the pulsars lie within the thick 
disk and their lines of sight do not probe structure behind 
them. Fortunately, this loss of resolution will not have 
a large influence on the results of this paper because we 
will find that the extragalactic contribution dominates 
the Galactic thick disk for £ > 5. 



B. Milky Way corona 

The Cordes & Lazio model does not include the Galac- 
tic corona component. Indeed, very few constraints on 
the electron density in the corona are available, which 
is unfortunate because the corona may contain a signif- 
icant fraction of the Milky Way's baryons j^. Models 
of the Milky Way suggest a baryonic mass in the bulge 
and disk components of ~ 5 x 1O 1O M0, whereas the virial 
mass is estimated at ~ 1O 12 M (2^. For the currently fa- 
vored baryon fraction fb = Qf,/Q m = 0.17±0.01 0, this 
suggests that an additional ~ 1.2 x 10 u Mq of "miss- 
ing" baryons should either be present inside the Milky 
Way's virial radius or else must have been ejected. (In- 
deed, Ref. j23| argued that at least half of the baryons 
are missing even though they assumed /& = 0.1. Increas- 
ing fb to 0.17 only makes the problem worse.) Assuming 
that a mass Mb of baryons were ionized and are present 
in a shell at radius R, they present a mean optical depth 



_ = (1 - Y/2)a T M b = 5 9 x iQ _ 5 (M fe /lO n M ) 



(.R/lOOkpc) 



(7) 



where Y — 0.24 is the Helium abundance, m p is the 
proton mass, and we have neglected the distance from 
the Earth to the center of the Galaxy in comparison with 
R. The virial radius for the Milky Way is estimated 
as r v i r ~ 260 kpc [23] , so that if a large fraction of the 
missing baryons are significantly inside the virial radius 
then the corona contribution to the optical depth could 



be comparable to the disk contribution. If instead the 
baryons are uniformly distributed in a ball, the radius 
and optical depth are related to the number density by 



R = 195 kpc 



t = 4.6 x 10" 



Mb/10 11 M Q 
?i H /10- 4 cm- 3 

M b \ 1/3 / 
l0 n M e ) V 



1/3 



"H 



io- 



2/3 

- 3 , •(«) 



One of the available constraints on the optical depth 
through the corona comes from the dispersion measures 
(DMs) of pulsars observed in the Large (LMC) and Small 
(SMC) Magellanic Clouds. The Thomson optical depth 
and the DM along a line of sight are both proportional 
to J n e dr, hence they are related by the equation r = 
2.05 x 10" 6 DM where the DM is in units of cm" 3 pc. The 
observed DM in the Magellanic Cloud pulsars is greater 
than that predicted by the Galactic thick disk model 0] 
by 16-91 cm -3 pc depending on the pulsar [2l]]. However 
much of this excess DM may be internal to the Magellanic 
Clouds and not part of the Galactic corona. Cordes & 
Lazio 21] argue for this interpretation since the DMs of 
pulsars not in the Magellanic Clouds are only fit if the 
electron density falls off beyond ~ 1 kpc from the plane of 
the Galaxy (see their Sec. 3.1.1). This is also consistent 
with the factor of ~ 5 disparity among DM excesses in 
the LMC along lines of sight separated by only a few 
degrees. In any case the contribution of the Galactic 
corona along the line of sight to the LMC at a distance 
of 50 kpc is limited to DM< 16cm~ 3 pc or r < 3 x 10~ 5 . 
If this upper limit is typical of lines of sight to 50 kpc and 
we use Eq. © with M b = 1.2 x lO n Af and R = 50 kpc 
as an upper limit on the optical depth beyond 50 kpc, 
we find that the optical depth through the Milky Way 
corona is at most 3 x 10~ 4 . If one instead assumes the 
uniform ball model (Eq. |SJ for the corona and assigns 
a baryonic mass M b = 1.2 x 10 n M Q , the DM to the 
LMC provides an upper limit of tt-h < 3.2 x 10~ 4 cm~ 3 
and r < 1.1 x 10~ 4 . This upper limit is significantly 
strengthened if the corona density is assumed to decrease 
outwards as a power-law more gradual than oc 1/r since 
this places the bulk of the baryons even farther away 
(if the power law is steeper than 1/r the optical depth 
diverges at small r, and would be cut off by either some 
core radius or our finite distance from the center of the 
Galaxy) . 

There are other constraints on the density of the Galac- 
tic corona. If the infall velocities of high-velocity clouds 
are assumed to be less than their terminal velocity due to 
drag in the Galactic corona, the density of the corona is 
limited to n H < 3 x 10" 4 cm" 3 at 100 kpc I2I Mod- 
els of the kinematics [25| and Ha emission|26] from 
the Magellanic Stream argue for nn ~ 10~ 4 cm" 3 and 
~ 5 x 10~ 5 cm~ 3 respectively at radius 50-65 kpc (but 
see criticism of these models arguing that survival of the 
stream then requires jih < 10 -5 cm -3 Interaction 
of high- velocity clouds with a Galactic corona is also sug- 
gested in order to produce some of the O VI detected in 
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Contamination of all-sky power spectrum 
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FIG. 1: The contributions to the B-mode polarization power 
spectrum. The white-noise (Ci tx I 2 ) lines are the lensing 
signal before cleaning (top line) and after cleaning assuming 
iterative lensing reconstruction with 0.25 fiK arcmin noise 
and 2 arcmin beam (bottom line). The dashed line shows 
the tensor B-mode spectrum for T/S = 10~°. The Galactic 
Thomson scattering contribution is shown by the solid line 
with points; note that due to the approximate parity sym- 
metry of the Galactic electron distribution, the odd I modes 
have much more power than even I. 



absorption [H H<|. If the higher values of n H [23. l29| 
are correct, and the corona remains this dense out to 
radii > 100 kpc, then it is possible that the optical depth 
through the halo may be ~ 10~ 4 . 

In order to assess the B-modes produced in the corona, 
we must remember that it is the fluctuations in optical 
depth St, rather than the mean optical depth f , that 
produces i?-modes. We will return to this point in the 
discussion, and for now we focus on the smooth Cordes 
& Lazio model for the known population of Galactic 
electrons. 



III. POWER SPECTRUM 

We next analyze the B-mode power spectrum of the 
scattered radiation. In the case of an all-sky experiment, 
one can simply do the spherical harmonic decomposition 



Cf B (scat) = 



1 



2£+l ^ 



J Y^(n)P(n)d 2 n 



(9) 



This is shown in Fig. ^ along with the lensing contri- 
bution, and the gravitational waves for T/S = 10~ 5 as 
computed by CMBFast |l7j . 

A realistic low-£ _B-mode experiment will exclude the 
Galactic Plane since this is where the foreground contam- 
ination is the worst. We would therefore like to under- 
stand the _B-modc contamination on the cut sky. Here 
we will consider cutting out the region within 10 degrees 
of the Galactic Plane; note that we are cutting this re- 
gion out of the final map P ga i > not out of the integration 



region in Eq. [5] Unlike the case of an all-sky map where 
Eq. © is the only "reasonable" (rotationally invariant) 
measurement of the _B-mode power spectrum, an analysis 
of the B-mode power spectrum on a cut sky requires the 
choice of an estimator. We choose the quadratic estima- 
tor, which involves the vector x = (Qi, ...Qn, Ui, ...Un) 
of length 2N, where N is the number of pixels. The 
covariance matrix C = (xx') is then written as 



Cij= N ij + g £ [(7f B Yf m (n i )YfJ(n i ) 



(io) 



where N is the noise covariance matrix, p a are the pa- 
rameters of the covariance matrix, and C Q are the power 
spectrum templates. Our objective is to estimate the 
parameters p a . The quadratic estimator determines the 
parameters using the relation 



Pa 

F al3 



I* \apqp, 

ixtC^C^x - ^Tr (Co^Cn^N) , 



Tr (Co x C a Co 'C/j) 



fill 



here F is the Fisher matrix and Co is a positive defi- 
nite Hermitian weighting matrix. The quadratic estima- 
tor method is unbiased regardless of the choice of Co in 
Eq. but is only optimal if Co is the true covariance 
matrix. (C controls the relative weighting of modes.) 
In our case, we will choose Co to contain the noise due 
to lensing after the iterative cleaning, so that it equals C 
if Njj is white noise with 0.83/xK arcmin and Cf B = 0, 
and Cf E >noise in Eq. (JTDJ since the CMB E-mode 
power spectrum is large. 

To understand the relationship between this method 
and the E/B decomposition of Ref. |30|, note that as 
Cf E — > oo, any polarization mode x that has non-zero 
inner product with an E mode has x^Cox — > oo. Also, 
any mode that is orthogonal to all E modes - i.e. that 
is a "pure i?-mode" in the terminology of Ref. |3(j - 
has Cox = Nx and so x is an eigenvector of Cq 1 with 
eigenvalue given by the inverse noise variance per pixel. 
Thus the quadratic estimator procedure Eq. (|llf> can be 
thought of as a method for projecting out the pure E- 
modes and ambiguous modes, and then computing the 
power spectrum; the Fisher matrix F serves also as the 
matrix of window functions for the quadratic estimators 
q a . (This same type of projection has been suggested by 
Ref. Hl|.) 

Numerically, we have evaluated Eq. {TJ using a weight 
with Cf E equal to 100 times the noise level for 2 < i < 
l max = 80. We use the implementation of the quadratic 
estimators described in Ref. |32^| (with no preconditioner 
for the Cq 1 operations) and the HEALPlX resolution 
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B-modes from Thomson scattering with Galactic Plane cut 
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FIG. 2: The B-mode power spectrum on a cut sky (with |6| < 
10 degrees removed) from Thomson scattering in the Galactic 
disk, computed with Eq. li lt . The dotted line is the cleaned 
lensing signal assuming iterative lensing reconstruction with 
0.25 /j,K arcmin noise and 2 arcmin beam, and the dashed line 
is the tensor contribution for T/S = 10~ 5 . The points show 
the polarization band powers from Thomson scattering (the 
open squares represent negative band powers). Some power 
may be missing at £ > 10 due to the sparse sampling of the 
pulsars used to construct the electron density model. 



(of order 1 kpc) is comparable to (and in some cases, 
greater than) the distance to the foreground-emitting re- 
gions. Since WMAP reports the temperature observed in 
a direction n but cannot specify the distance to the emit- 
ting region, one needs an additional assumption in order 
to calculate P ga i ■ We will make the assumption here that 
the microwave foreground can be decomposed into a "lo- 
cal" contribution from nearby (distances much less than 
the size of the Galaxy) and a "global" contribution from 
distances of order the size of the Galaxy. The local con- 
tribution we assume to be plane-parallel; in this case the 
quadrupole is independent of the position of the scatter- 
ing electron and we may use the observed quadrupole at 
Earth. The global contribution to the photon quadrupole 
at the position of the scattering electron can also be taken 
to equal the foreground quadrupole observed at Earth, 
since the distance to the scattering electron (< 1 kpc) 
is small compared to the size of the Galaxy. Being non- 
axisymmetric, the global contribution can produce both 
E and B polarization. Within this set of assumptions, we 
can calculate the i?-modes from Thomson scattering of 
Galactic foreground emission by repeating the procedure 
of Sec. IIIII with the uncleaned WMAP maps. Since these 
maps contain both foregrounds and CMB, the result - 
shown in Fig. [3] - is the B-mode power spectrum from 
Thomson scattering of both the CMB and the Galactic 
foreground emission. 



5 pixelization system [3i| , which has 12 288 pixels. In 
Fig.|21we have shown the power spectrum of the B-mode 
polarization obtained if we accept only the 10 112 pixels 
at Galactic latitude |6| > 10 degrees. We have plugged 
the resulting polarization data vector x of length 20 224 
into Eq. JTT}, leaving out the Tr (Cq 1 C a C^" 1 N) ter m m 
q a since the Thomson scattering map does not include 
the noise. The resulting C^f B (scat) represent the addi- 
tive bias to the -B-mode power spectrum estimators due 
to the Galactic Thomson scattering component. (These 
are negative for £ = 4, 5; note that on a cut sky it is possi- 
ble for the power spectrum estimator to become negative 
even without noise subtraction.) 



V. EXTRAGALACTIC CONTRIBUTION 

In addition to the Galactic electrons, the CMB 
quadrupole can also be scattered off of nearby struc- 
tures in the universe, generating B-mode polarization. 
In this section, we begin by computing the B polariza- 
tion from structures out to a distance of 80 hr 1 Mpc 
using a model for the local electron distribution. We 
then calculate the expected Cf B statistically from the 
power spectrum. The former method has the advantage 
of corresponding approximately to the actual realization 
of the electron distribution in the universe, whereas the 
latter method is better for addressing the contribution 
from larger distances. 



IV. POLARIZATION FROM SCATTERED 
GALACTIC EMISSION 

In addition to scattering the CMB quadrupole, Galac- 
tic electrons can scatter the microwave foreground emis- 
sion emitted by the Galaxy itself. In order to compute 
this effect, one needs not only a model for the distribution 
of electrons but also for the foreground emission, so that 
one can calculate the quadrupole in Eq. JBJ. While the 
foreground emission at frequencies from 23-94 GHz has 
been mapped out by WMAP, the quadrupole T 2m (rn) 
at the scattering electron may differ from the locally ob- 
served quadrupole because the distance to the electron 



A. Constrained realization method 

Our first method for estimating the local extragalac- 
tic contribution to the B-mode polarization is to use 
the "constrained realization" A^-body simulations from 
Mathis et al. |3J| and assume that the electron den- 
sity traces the dark matter. These simulations have 
constrained initial conditions intended to reproduce the 
structures observed in the local universe out to a distance 
of ~ 80/i _1 Mpc and on scales larger than 5ft. -1 Mpc (on 
smaller scales, the initial conditions are filled in with a 
Gaussian random field with appropriate power). Mathis 
et al. [34| simulated both a ACDM and Einstein-de Sitter 
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(a) B-modes from Thomson scattering in K-band 
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(b) B-modes from Thomson scattering in W-band 
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FIG. 3: The B-mode power spectrum on a cut sky (with \b\ < 10 degrees removed) from Thomson scattering, in (a) the K 
(23 GHz) band and (b) the W (94 GHz) band. This figure is the same as Fig. [5] except that we have replaced the cosmic 
quadrupole with the raw K and W band quadrupoles observed by WMAP (no frequency cleaning) . Some power may be missing 
at £ > 10 due to the sparse sampling of the pulsars used to construct the electron density model. 



cosmology; here we have used the former (with fi ro = 0.3, 
Ha = 70 km/s/Mpc, and a% = 0.9) as it corresponds 
most closely with the currently favored cosmological pa- 
rameters. 

The perturbation to the optical depth due to nearby 
extragalactic structures is calculated according to 



wavenumber, 



5t 



ax 



n e (z) 



1 



8 e dr, 



(12) 



where 5 e is the fractional electron density perturbation; 
here we assume S e = 6, which should be valid except on 
very small scales where the baryons are segregated from 
the dark matter. We have written this as an integral 
over the comoving distance r. The mean electron density 
assuming complete ionization is 



_ 3Sl b H$(l-Y/2) n 



2.2xl0~ 7 (l + z) 3 cnr 



(13) 

where G is Newton's gravitational constant and the 
baryon density Vl^h 2 is from WMAP Q. Here we have 
neglected redshift evolution of n e (z)/(l + z) and S e since 
the simulation goes out to a distance of 80/i _1 Mpc 
(z = 0.027); n e (z)/(\ + z) is only 5.4% greater at 
z = 0.027 than today, while S e is slightly less (by 1.3% in 
the linear regime) . Equation i|12[l cannot be used directly 
because the dark matter distribution in the simulation 
is represented by individual particles; we have therefore 
smoothed with a Gaussian of la width 0.1/i -1 Mpc in 
each dimension. The Poisson noise at the smoothing scale 
k = 10/i/Mpc is 4% of the linear matter power spectrum. 

Physically, the free electrons cannot trace the dark 
matter down to arbitrarily small scales, and one would 
expect smoothing on scales shorter than the Jeans 



kj= v L = 4,2ft Mpc 1 J——(1 + S), 

C S V Tg 

(14) 

where c s is the sound speed, T g is the gas tempera- 
ture, and we have assumed the sound speed relation 
c 2 s = 5/csT/3/i with mean mass per particle [i = 0.59m p , 
appropriate for ionized gas with Y = 0.24. The inter- 
galactic medium is photoionized so we expect T g > 10 4 K 
and hence in regions of average density all structures 
smaller than kj 1 > 0.2/i _1 Mpc should be washed out if 
we consider 5 e instead of S. In this sense our estimate us- 
ing smoothing at O.lh^ 1 Mpc is conservative and should 
overestimate Cf B . Most of the small-scale power is how- 
ever in overdense regions where 1 + <5 ^> 1 , so depending 
on T g in these regions the smoothing scale may be dif- 
ferent from the typical value of fcj 1 . We find that the 
£ < 10 modes are only slightly affected (by < 20% ) if we 
increase the smoothing length to O.bh' 1 Mpc, indicat- 
ing that the results are not dominated by the small-scale 
modes. 

Once St maps are constructed, these can be fed 
through Eq. © to compute the polarization signal. The 
power spectrum on a cut sky (retaining the |6| > 10 de- 
grees region) can then be determined from a quadratic 
estimator, as was done in Sec. IIIII 



B. Power spectrum method 

Our second approach to computing Cf B is to estimate 
the power spectrum of the electron density perturbations, 
and then convert this via Eq. JSJ into a power spectrum 
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for the polarization. This approach is purely statistical 
and does not require any particular realization of the 
electron distribution. Its principal disadvantage is that 
the sky r(n) may look very different from different parts 
of the universe, and thus the local value of CJ T may 
deviate substantially from the global average computed 
here. 

The perturbation to the Thomson optical depth in di- 
rection n is given by Eq. (|12|l . We begin with the case 
where 5 e is a Fourier wave, S e (r, n) = e"" kn ~ er , where e is 
an infinitesimal parameter that forces convergence of the 
integral at infinity. We then have i5r(n) = n e <7T / '{e—ikji) 
where fi = k • n is the angle between the wavevector and 
the line of sight. Taking the spherical harmonic trans- 
form, we find 



o 

iSmO 

2k 



d(f> J d^i F fm (n)r(n) 
n e a T V^WTT) f (15) 



where Pi is a Legendre polynomial and we have taken k 
to be along the z-axis. If we define the coefficients 



at 



1 Pe(n)dfi 



then using the recursion relation 



(16) 



tP t (p) = {2£- l)nP t -i{l*) ~ {l-l)Pt-2{fi), (17) 

and noting that for I > orthogonality implies 
j^Pi^fijd/i = 0, we find ap = — (£ — l)ai~2/£- It is 
then easy to show by induction from the initial values 
clq = —iir and ai = 2 that 



2-V- 1 7rr(£+ 1) 



1 



(18) 



To generalize from the case of a single Fourier mode 
to a random field, we need to incoherently integrate the 
power spectrum CJ T over all Fourier modes: 



cr 



21+1 

7T 3 [r(l+l)] 2 (fMT T ) 2 

4'{r[(*/2) + 1]} 4 



dk 
~~k 



dk 



k-'Ai(k) T ; (19) 



in the last line we have substituted the result of Eq. (|18fl . 
(The summation over m is rotationally invariant and so 
it does not matter that we have assumed k to be in the 
z-direction.) This is slightly smaller than the Limber 
approximation 



II T 



Limbc 



f _ 27T 2 

(n e a T f / r 2 -^-A 2 Sc (k)\ k=(e+1/2)/r dr 
2ir 2 (n e a T ) 2 



£+1/2 



(20) 



by only 5% for £ = 1 and 2% for 1 = 2. 

To go from Eq. (|19|l to the B-mode power spectrum, we 
note that from Eq. J^J), the multipoles of the polarization 
from the local universe are 



B, 



loc.im 



|Yfi(n)P(n)d 2 n 

Tit." rp ^ 

/ , 1 mm'm" - L 2m''l"m" t 
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(21) 



where we have used the mode coupling integral 

Cl>m» = J YfJ (n) Yf m , (n)Y e „ m „ (n) d 2 n. (22) 

Expressing the tensor spherical harmonics in terms of the 
spin- weighted spherical harmonics, 

[Yf m ] q = [Yf m ] u = \{Y 2 n + Y e - 2 ) and 



[ Y £m] U 



[ Y £m] Q ~ Yi ( Y t™ ~ Y im)' 1 - ?> > 



gives 



Til" _ _L / f V -2* v 2 _ v 2* v -2\ v ,2- 

1 mm'm" ~ r>„- / X 1 £m 1 2m' 1 £m 1 2m') * m " 

> h{2i + 1){2£" + 1) 



2i 



-i(-iy 



1 l" 2 

2 0-2 



An 

£ £" 2 



n / / (24) 
mm m ' 



for I" = £ + 1 and l" m , m i, = otherwise by symmetry 
(for £" — t = or ±2, I^ m , m „ = vanishes by parity, for 
\i" — £\ > 2 the 3j symbols vanish). Here we have used 
the three spherical harmonic integral 0, l36l . l3?i| . Now 
the power spectrum of the Thomson-scattered radiation 
is obtained by taking the mean square value of Eq. 121fl : 

Cf B (l0c) = 1 £ \Cm'm"\ 2 (\T2m'r e "m"\ 2 ) 
m' £" m" 

= ~^Q^2 T l^mm'm"| 2 CJ"' (25) 

m' £"m" 

where we may drop cross-terms between different values 
of m' , £", and m" since the different multipoles are uncor- 
rected and we have assumed that r in the local universe 
is independent of the CMB quadrupole. (The value of 
the sum cannot depend on m because of rotational sym- 
metry.) We may do the summation over ml and m" by 
applying the 3j symbol orthogonality relations, yielding 
(for t" = t + 1) 

5(21" + 1) f £ £" 2 x 2 
4tt { 2 -2 



V \i u " ' A 

/ j I mm m I 



(26) 



Using the tabulated values of this specific form of the 3j 
symbol [36j gives the final result 



Cf B (loc) = 



iCi 1 [{l + 2)CV_ 1 + {£-l)Cr +1 ] 
80tt(2£+ 1) 



(27) 
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Results 



statistically using its covariance matrix 



We can evaluate Eq. l(2"Tj> given any power spectrum 
A 2 (fc), which here we will assume to be the nonlinear 
matter power spectrum. We have computed the spec- 
trum using the Eisenstein & Hu [3^| transfer function 
and the Peacock & Dodds nonlinear mapping. The 
results from this calculation are shown in the thick line in 
Fig-EIa). For comparison, we have also shown (thin line) 
the results obtained by cutting off the integral in Eq. l(T^|l 
at fc min = (£+l/2)/r max , where r max = 80/i _1 Mpc is the 
radius to which we integrate in the constrained realiza- 
tion method. Within the context of the Limber approx- 
imation, this measures the amount of this power that is 
recovered by the constrained realization method. The 
difference between the thick and thin curves in Fig. 0Ja) 
represents the contribution to Cf B coming from struc- 
tures farther away than 80/i _1 Mpc; this is 50% of the 
power at t = 10, with nearby structures dominating at 
lower £ and more distant structures at higher I. 

Also shown is the power spectrum from the constrained 
realization method (points in Fig. ^i). We can see from 
the figure that the power spectrum Cf B obtained from 
this method is greater than the predicted result out to 
a distance of ^Qh" 1 Mpc (thin line). This indicates that 
the nearby universe is more inhomogcneous than aver- 
age, which is not entirely surprising since we live in an 
overdense region of the universe and it is these regions 
that contribute most of the small-scale power. 

Since the Mathis et al. |3J| simulations provide an 
actual optical depth map, it is possible to construct an 
actual realization for the combined Thomson scattering 
signal from the Milky Way's disk and from nearby ex- 
tragalactic structures. We have computed the power 
spectrum Cf B of this realization with the quadratic 
estimator, and added to it the "missing" power from 
r > 80/i _1 Mpc; the result is shown in Fig. Efb). This is 
our final result for the frequency-independent, Thomson- 
induced £?-mode power spectrum from the Galactic disk 
and nearby extragalactic structures (but excluding the 
corona). The contribution is roughly (1-2) x 10 _4 /xK per 
logarithmic interval in i, which is well below the cleaned 
lensing signal even for an optimistic experiment. 

We can be more quantitative about the effect on the 
estimation of the tensor signal T/S by using the ten- 
sor signal C T /5 = dC/d(T/S) as a template in the 
quadratic estimator instead of the band powers. In 
this case we can get the spurious contribution to T/S 
from the equation A(T/5) = (It/s/Ft/S ,tJS i where qr/s 
and Ft/s,t/s are computed using Eq. i|ll[) with Cjyg. 
The map x to be inserted into the equation for qr/s is 
the sum of the known Galactic disk and extragalactic 
(r < SOft^Mpc) signals and the unknown extragalac- 
tic signal from r > 80/i _1 Mpc, which we denote 
and x^ 2 ) respectively. We can put into the estimator the 
actual value of x^ 1 ), whereas we must account for x' 2 ) 



A(T/S) = 



(Qt/s) 
Ft/s,t/s 

(Ft/s,t/s) 1 
1 



Tr(C^C T/s C^C( 2 ); 



(28) 



This equation yields A(T/S) = 1.0 x 10~ 7 for the tem- 
plate with optical depth r = 0.17, and 3.6 x 10~ 7 for 
t = 0.07. This is negligible, since even for the optimistic 
(0.25 /iK arcmin noise, 2 arcmin beam) experiment, the 
Fisher matrix uncertainties in T/S are 6.2 x 10~ 7 and 
1.7 x 10~ 6 for t = 0.17 and 0.07 respectively. 



VI. DISCUSSION 

We have considered the contribution to the CMB B- 
mode polarization from Thomson scattering in the local 
universe, which is a potential contaminant to the gravi- 
tational wave signal. Our findings are that: 

1. The Thomson scattering signal from the Galactic 
disk and from nearby extragalactic structures is 
well below the weak lensing B-modes even after 
the lensing signal is "cleaned" out using a lensing 
reconstruction from a very optimistic experiment. 
This signal is believed to be negligible. 

2. The Galactic disk also Thomson-scatters the fore- 
ground radiation emitted by the Milky Way. At 
the low multipoles, this radiation can exceed the 
cleaned lensing signal; in the K band it is roughly 
equal to the uncleaned lensing signal at I = 2- 
3. However this signal is strongly frequency- 
dependent, and in principle can be removed by the 
same techniques used to clean other Galactic fore- 
grounds. 

3. Thomson scattering within the hot Galactic corona 
makes an unknown contribution to the i?-modes, 
depending on the flucutations in the optical depth 
through the corona. The contribution can be 
roughly estimated based on Eq. (|27|l : 



1(1+1) 
2tt 



a 



BB 



80tt t 



2.4 iiK 2 A 2 .; 



(29) 



this can exceed the cleaned lensing signal if the 
large angular scale fluctuations in the optical depth 
through the corona are A T > 2 x 10~ 4 . Since a 
mean optical depth through the corona of up to 
f ~ O(10~ 4 ) appears to be reasonable, and there 
are only very limited constraints on the angular 
distribution of free electrons in the corona, fluctu- 
ations of A T > 2 x 10~ 4 cannot yet be ruled out. 
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(a) B-modes from nearby extragalactic structure (b) Total (Galactic + extragalactic) B-modes 




2 3 4 5 10 15 20 30 2 3 4 5 10 15 20 30 

Multipole, L Multipole, L 

FIG. 4: The contribution from nearby extragalactic structures to the .B-modes. Panel (a) shows the power spectrum for 
the constrained realization method of Sec. IV Al (points') . the power spectrum method of Sec. IV Bl (thick line), and the power 
spectrum method with the integration cut off at r max = 80/i _1 Mpc (thin line). Panel (b) shows the total B-mode from 
local scattering (points), constructed by taking the optical depth map from the Galaxy and nearby structures, computing the 
resulting Cf B using the quadratic estimator, and then adding the estimated power from r > SOh 1 Mpc. For reference, the 
cleaned-lensing and tensor curves from Fig. H are shown. 



Of these, the scattering from the Galactic corona is of 
the most interest. In principle, detection of this signal 
would provide evidence that the corona contains some of 
the missing baryons; non-detection of the B-mode sig- 
nal would provide an upper limit on the anisotropy of 
the corona. However the practical difficulties should not 
be understated: reaching the "cleaned lensing" limit of 
Cf B = 5.8 x 10 _8 /iK 2 requires a detector sensitivity of 
0.25 /iK arcmin, two orders of magnitude better than the 
upcoming Planck satellite 0] • It * s a l so possible that the 
frequency structure of the polarized foregrounds is suffi- 
ciently complicated that extraction of small large-angle 
B-mode signals becomes impractical. Finally it is pos- 
sible that tensor perturbations, or some other source of 
frequency-independent i?-modes such as cosmic strings, 
are present and swamp the Galactic Thomson scattering 
signal. 

From the perspective of CMB gravitational wave 
searches, the Galactic Thomson scattering signal appears 
to be negligible for all except the most futuristic exper- 
iments. Part of the reason for this is the low value of 
the local quadrupole Cj T , which is only ~ 20% of the 
ACDM expected value; if ACDM is indeed correct, we are 
in this respect simply lucky. For the very futuristic ex- 



periments, however, if a frequency-independent £?-mode 
signal were observed close to the detectability threshold, 
then the interpretation as a gravitational wave source 
would require that Thomson scattering in the local uni- 
verse be robustly ruled out as the source. This would 
require significant improvement in our understanding of 
the Galactic corona. 
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